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ABSTRACT 


We use a Green’s function technique to calculate the position of deep defects in 
narrow gap semiconductors. We consider substitutional (including antisite), vacancy, 
and interstitial (self and foreign) deep defects. We also use perturbation theory to look 
at the effect of non-parabolic bands on shallow defect energies and find non- 
parabolicity can increase the binding by 10% or so. We consider mercury cadmium 
telluride (MCT), mercury zinc telluride (MZT) and mercury zinc selenide (MZS). For 
substitutional and interstitial defects we look at the situation with and without 
relaxation. 

For substitutional impurities in MCT, MZT, and MZS, we consider x (the concentration 
of Cd or Zn) in the range 0.1 < x < 0.3 and also consider appropriate x so E g = 0.1 eV 

for each of the three compounds. We consider several cation site s-like deep levels 
and anion site p-like levels. For E g = 0.1 eV, we also consider the effects of relaxation. 
Similar comments apply to the interstitial deep levels whereas no relaxation is 
considered for the ideal vacancy model. Relaxation effects can be greater for the 
interstitial than the substitutional cases. 

Specific results are given in figures and tables and comparison to experiment is made 
in a limited number of cases. We find, for example, that I, Se, S, Rn, and N are 
possible cation site, s-like deep levels in MCT and Zn and Mg are for anion site, p-like 
levels; both levels for substitutional cases. The corresponding cation and anion site 
levels for interstitial deep defects are (Au, Ag, Hg, Cd, Cu, Zn) and (N, Ar, O, F). For the 
substitutional cases we have some examples of relaxation moving the levels into the 
band gap, whereas for the interstitial case we have examples where relaxation moves 
it out of the band gap. Future work involves calculating the effects of charge state 
interaction and seeing the effect of relaxation on vacancy levels. 
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L INTRODUCTION 


In the semi-annual report* 31 ), we have already presented the background for this work 
so this report will mainly focus on results. 

The main defects we will look at are deep defects caused by short range potentials. 
We will include substitutional (including antisite), vacancy, and interstitial (self and 
foreign) deep defects. Some brief results will be given for shallow defects caused by 
long range potentials, where perturbation theory will be used to look at the effects of 
non-parabolic bands. 

We will consider mercury cadmium telluride, mercury zinc telluride and mercury zinc 
selenide. They have non-parabolic bands, but are interesting because they have 
large mobilities, tunable band gaps, are candidates for microgravity growth and are 
potentially important materials to make infrared detectors* 32 ' 34 ). 

As is standard, the Green's function technique will be used to calculate the positions of 
the deep defects. 

Note that deep defects may or may not be near the band edge. They are often very 
important as recombination centers. Shallow defects are near the band edge and in 
narrow gap semiconductors are often assumed to be fully ionized. Their importance is 
to supply carriers, but like all defects they may also scatter carriers. 

II. DEEP LEVELS OF POINT DEFECTS IN SEMICONDUCTORS 

A. Green’s Function Method 


For substitutional defects we showed in the semi-annual report* 31 ) that the 
Schrodinger wave equation and the definition of Green’s Functions lead to 


det 


li-Gi (E)V 


= 0 


(D 


which determines the defect energy levels. This equation is true for both 
A-i (s-like) and T 2 (p-like) states, with appropriate choice of V { and G° (E) as 
discussed in the semi-annual report. The relevant Green’s Functions are 
Gj (E) and the potential of the defect determines the V, 
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It turns out that the energy bands can be written in terms of diagonal matrix 
elements as 


levels: G aa = Ug 
T 2 levels: G aa = Up 1 

where U s and U p are impurity potentials determined by V s and V p . 

B. Ideal Vacancy Model 


(2a) 

(2b) 


For the ideal vacancy mode! U s and U p can be assumed infinite which 
equivalently cuts all the bonds, and so 


Gaa( E ) = 0 

Gaa(E)= 0 

C. Band Structure 


(3a) 

(3b) 


In order to do our calculations we need to know the band structure, 
because G° , the Green’s functions of the perfect crystal, depends on the 
band structure. For ll-VI semiconductors, in principle we also need to include 
spin orbit interactions. When they can be ignored, which may happen, the 
presence of spin just doubles the degeneracy. 


D. Lattice Relaxation 


Lattice relaxation, particularly for the substitutional case, has been discussed in 
the semi-annual report. We showed that differing host and impurity body 
lengths, due to relaxation, determined the off diagonal elements of the defect 
potential. The amount of relaxation was computable by molecular dynamics. 
The ideal vacancy mode! does not allow relaxation, and the peculiarities of 
relaxation due to interstitials will be discussed in the next section. 


2 



E. Interstitial Defects 

Interstitial impurities play an important role, as point defects, in semiconductors. 
In fact, self- interstitial impurities and vacancies are the most important native 
defects in semiconductors. From experiment, it is found that Li, Mg, and the 
transition metals in Si prefer to occupy the interstitial sites. Furthermore, in the H- 
VI materials, it is believed that substitutional group-1 acceptors may 
spontaneously move to an interstitial site< b1 ). In this section, we study tetrahetral 
- site interstitial s - and p - bonded impurities in MCT, MZT, and MZS. We 
believe that hexagonal interstitial impurities can also be calculated using a 
similar model to that described below. 

For a substitutional defect in zincblende structure semiconductors, the defect 
has a tetrahedral symmetry in the materials. That is, there are four nearest- 
neighbor atoms and twelve second-nearest-neighbor atoms around the defect. 

If the distance between the defect and the nearest neighbors is d, then the 
distance between the defect and its second neighbors is 1 .64d. However, the 
situation of tetrahedral-site interstitial impurities is quite different. In this case, 
there are four nearest neighbors around an interstitial impurity with a distance d, 
but six second neighbors around it with an only slightly greater distance 1.1 5d. 
Hence, the nearest-neighbor approximation, which is often used for 
substitutional defects, is not acceptable for interstitial impurities. In contrast, in a 
crude approximation the interstitial may be thought of as tenfold coordinated. 
That is, we may assume that there are ten nearest neighbors around the 
interstitial impurity. 

Sankey and Dow( b1 ) proposed a model to calculate deep levels introduced by a 
s- and p-bonded impurity occupying an interstitial site of tetrahedral symmetry in 
Si. This model can be easily extended for ll-VI materials. The main points of this 
model are as follows: The host Hamiltonian H 0 is determined by the Vogl’s 

ION x ION sp 3 s band structure theory, where N is the number of unit cells, and 

the ten states are sT, si , pxt, pxi, pyT, pyi, pzT, pzi, s*T, and s*i . The use of 
excited s* orbitals allows a description of the conduction bands that is flexible 
enough to treat both direct and indirect band gap semiconductors. The 
interstitial impurity is treated by the ION host basis orbitals and includes one s 
orbital and three p orbital centered at the interstitial site. Hence the Hamiltonian 
matrix for the host plus interstitial impurity is a direct sum H 0 and Hj, where H ( is 
the diagonal 4x4 interstitial Hamiltonian. 
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Therefore, the modeled Hamiltonian takes the schematic block form 


H = 


H 0 W 

wV 


( 4 ) 


where W denotes the coupling matrix elements between the interstitial atom 
and the host. 


The calculation of H 0 is exactly same as that in the case of substitutional 
defects. H 1 is a 4 x 4 diagonal matrix which has a form as 



(5) 


where e s and e p are the orbital energies of the interstitial impurities. For 

example, for a neutral self-interstitial impurity in Si, e s and e p are the orbital 

energies of Si. The elements of the transfer matrix W are taken to be zero 
except those between the interstitial atom and its first and second neighbors. 
The elements coupling the interstitial and the nearest neighbors, W nn1 , are 

taken to be those of the Vogl et al.( b2 ) tight - binding model and the elements 
coupling the interstitial and the second neighbors, W nn2 , are taken to be a 

fraction of W nn1 . In the case of Si, Sankey and Dow assumed that 

^nn2 = O-^Wnnl (6) 

The deep levels are calculated using the Green’s function method as usual. 
Therefore, they are determined by the determinanthal equation 


det[i - G ref (E) V] = 0, (7 ) 

where G ref (E) = (E - H ref ) _1 , and H ref is the Hamiltonian of eq (4) with H 1 = H' ref . 
Here H ! ref is the Hamiltonian of the neutral, self-interstitial impurity. V is the 
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defect potential matrix which has the same form and meaning as that in 
Hjalmarson’s theory. Since the interstitial atom occupies a site of tetrahedral 
symmetry, the point group symmetry is T d , The condition for a deep level in the 
band gap simply reduces to the familiar conditions as below: 

[G r e f (E)] ss =(U s )- 1 1 

[Gref(E)]pp = (Up)' 1 . (8) 

To extend this model to be used for ll-Vi materials, we just need to include the 
spin-orbital coupling in the calculations. In the new model, the host Hamiltonian 

H 0 is now determined by the 20N x 20N sp 3 s* (spin with spin orbit coupling 

doubles the number of states) band structure theory proposed by Kobayashi( b3 ), 
and the rest of the model will be the same as Sankey’s model. 

F. Summary of Results for MCT, MZT, and MZS 

1. Results for Substitutional Impurities (including antisite defects) 

(a) Deep levels as a function of the composition x 

The results for the chemical trends of deep levels produced 
by various impurities in MCT, MZT, and MZS as a function of 
the alloy composition x are summarized in Fig. 1 - Fig. 5. 

Fig. 1 shows the cation - site, s - like deep levels produced 
by N, S, and I impurities in MCT as the function of x in the 
range of 0.1 to 0.3. Fig. 2 shows the anion - site, p - like 
deep levels produced by Mg, Zn, and Cd impurities in MCT 
for x in the same range. Fig, 3 to Fig. 5 show deep levels in 
MZT and MZS for x from 0.1 to 0.3. It is worth mentioning 
that there are no cation-site, s-like deep levels predicted in 
MZS. 

One of the most important conclusions from these figures is 
that the slopes of the deep levels dE/dx are always smaller 
than dEgap / dx. This usually has the consequence that, if a 
deep level exists for x = 1 , then this deep level will remain in 
the band gap for a very large range of x. It is also worthy to 
note that the slopes of the anion - site deep levels are very 
small, these levels, thus, are “attached” to the valence band. 
These two results are very useful in practice. For example, 
detecting impurity energy levels in narrow gap alloys is 
much more difficult than that in wide gap alloys. However, 
following the results above, instead of detecting the deep 
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levels in the narrow gap alloy, we can try to detect the same 
deep levels in the same alloy with a wider gap. If we find an 
anion - site deep level in a wide gap alloy, we can be 
almost sure that this level will be found in the same alloy 
with a narrow gap. 

(b) Predicted deep levels with no relaxation 

The results of the predicted deep levels produced by 
various impurities in MCT, MZT, and MZS are summarized 
in Fig. 6 to Fig. 1 0. The effect of lattice relaxation on deep 
levels is not included in the present predictions but will be 
included later. The energy band gap of the materials is 
specified as 0.1 eV. Thus, the compositions of MCT, MZT, 
and MZS considered here are approximately 0.22, 0.15, 
and 0.08 respectively, which will give the desired band gap 
at T = 0 K (approximately). Fig. 6 and Fig, 8 show the cation 
-site, s - like deep levels, and Fig. 7, Fig. 9, and Fig. 10 show 
the anion - site, p - like deep levels. In these figures, 
besides the deep levels within the band gap we also show 
the deep levels resonant near the band gap edges. 

Inspecting these figures, we can see that all impurities 
substituted for cation atoms (Hg, Cd, and Zn) in MCT (Fig. 6) 
and MZT (Fig. 8) are donor - like impurities. That is, they are 
all from the rows at the right side of the cation atoms in the 
Periodic Table. Therefore they are more electro- negative 
than the cation atoms in MCT and MZT. On the other hand, 
all impurities substitutional for anion atoms (Te and Se) of 
MCT (Fig.7), MZT (Fig. 9), and MZS (Fig. 10) are acceptor 
- like impurities. They are all from the rows at the left side in 
the Periodic Table, and they are, thus, more electropositive 
than the anion atoms in MCT, MZT, and MZS. This result 
also has its practical use. Typically, if we find an impurity 
level produced by substitution for cation atoms in MCT and 
MZT, then this impurity level is donor - like. On the other 
hand, if we find an impurity level produced by substitution 
for anion atoms, then this level is acceptor - like. 

It can be also seen from these figures that even though the 
same impurities will form different energy levels in different 
materials, the order of these levels in the materials is the 
same. For example, the order of deep levels produced by 
impurities of Mg, Zn, and Cd substituted for anion atom Te 


6 



in MCT is E Mg < E Zn < E Cd . The same order is found in 

MZT and MZS. Therefore if we can determine one deep 
level produced by a specific impurity from experiments, we 
can then guess the positions of deep levels produced by 
other impurities. This illustrates what we mean when we 
say our calculation is very good for predicting chemical 
trends. 

(c) Effects of lattice relaxation 

The results of the deep levels in the band gap as functions 
of the appropriate on - site defect potential, U A1 or U T2 (see 
section IIA ), denoted as “impurity potential”, and the 
appropriate off - diagonal parts of the defect potentials, a A1 
or ( see al), are summarized in Fig. 1 1 to Fig. 15 
in MCT, MZT, and MZS. The curves in these figures 

correspond to the various values of a shown in the figures. 
As mentioned before, the effects of lattice relaxation are 

included in parameter a, and a = 0 corresponds to the 
case of no lattice relaxation. In the case of MCT in Fig. 1 1, 
where we consider the cation - site, s - like ( A-j ) deep 

levels, a A1 = 0.5 corresponds a 27% outward relaxation, 

and a A1 = -0.5 corresponds a 15% inward relaxation. In 
Fig. 12, where we consider the anion - site, p -like ( T 2 ) 

deep levels, a T2 = 0.2 corresponds a 12% inward 
relaxation and a T2 = -0.2 corresponds a 1 8% outward 

relaxation. It is worthy to note that the relaxation directions 
are opposite for A 1 and T 2 states. Similar calculations can 
be done for MZT and MZS in Fig. 13 to Fig. 15. 

The deep levels for a particular impurity with a particular 
value of a can be obtained from these figures by finding the 

intersections with the curve for that value of a with a vertical 
line drawn from the value of “ impurity potential “ for that 
impurity. If there is no intersection for a particular case, 
then no deep level of the appropriate symmetry is predicted 

in the band gap for the impurity with that value of a. 

Inspection of these figures reveals the qualitative effects of 
lattice relaxation on deep levels. First, these effects are very 
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small for cation - site, s - like deep levels in MCT and MZT. 
They can be totally neglected in the calculations of deep 

levels for specific impurities. It can be also seen that the a * 
0 curves all have chemical trends which are similar to the 

a = 0 ( no lattice relaxation ) curves. From Fig. 12 we can 
see that an inward relaxation (a-j- 2 > 0 ) moves deep levels ' 

toward the valence band and an outward relaxation ( a T2 < 

0 ) moves deep levels toward the conduction band. Since 
the band gap we consider here is only 0.1 eV, the effects of 
lattice relaxation on deep levels are very important. 

Including these effects may move deep levels out of the 
band gap or move deep levels into the band gap. For 
example, in Fig. 14, an impurity with “impurity potential” = 

2.9 eV (Mg) has a deep level in the middle of the band 

gap for a = 0 (no relaxation ). However, a 10% inward 
relaxation will move this energy level out of the band gap, 
the energy level will resonate in the valence band. On the 
other hand, a 15% outward relaxation will move this energy 
level to the conduction band. 

(d) Effects of lattice relaxation on deep levels 

In order to determine a deep level for a particular impurity in 
MCT, MZT, and MZS with appropriate symmetry, a must be 
specified. To specify a d ( and d H must be determined first. 

The host bond length d H as a function of the composition in 

MCT, MZT, and MZS can be calculated using d H = r^+ r^ 

> 

and the impurity bond length d|Can be determined by 

molecular dynamics technique discussed in section IID. The 
effects of lattice relaxation on deep levels produced by 
various impurities in MCT, MZT, and MZS are summarized 
in Table 1. Since these effects on cation-site, s-like levels 
can be ignored, only the results on anion-site, p-like 
levels are shown in Table 1. 

It can be seen from Table 1 that the percentage of lattice 
relaxation is in the range of 4% - 6.5%. The shift of deep 
levels after considering the relaxation is less than Eg/2 or 
0.05 eV. Therefore, if a deep level with no relaxation is in 
the middle of the gap, this level will stay in the gap after 
relaxing. However, if a deep level with no relaxation is near 
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the gap edges, this level may move out of the gap after 
relaxing, or move in the gap from the conduction band or 
the valence band. For example, antisite defect Zn in Te site 
in MCT is predicted to produce a deep level in the middle of 
the band gap (Eg/2) without relaxation, and including the 
effects this level is predicted to move to near the conduction 
band (0.9 Eg). When Tl and In are substituded for Se in 
MZS, they are predicted to form deep levels just below the 
valence band without relaxation, but including the effects, 
these two levels move into the band gap. 

2. Results for Vacancies 

The deep levels produced by vacancies in MCT, MZT, and MZS are 
summarized in Fig. 16 and Fig. 17. Fig. 16 shows the cation-site, s-like 
vacancies levels of MCT, MZT, and MZS. The vacancy levels of MCT and 
MZT are resonant in the valence band and the vacancy level of MZS is in 
the middle of the band gap. Fig. 17 shows the anion-site, p-like vacancy 
levels. These levels are slightly below the corresponding levels in Fig. 

16. These levels are calculated by “the ideal vacancy model” mentioned 
in Sec. IIB. Since the effects of lattice relaxation are not included in the 
calculation, the accuracy of these levels need to be improved. This 
problem will be discussed in the last chapter. 

3. Results for Interstitial Impurities 

(a) Predicted deep levels with no relaxation 

The results of the predicted deep levels produced by interstitial 
impurities in MCT, MZT, and MZS are summarized in Fig, 18 to Fig. 
23, Fig. 18, Fig.20, and Fig. 22 show the cation-site, s-like deep 
levels, and Fig. 19, Fig, 21, and Fig.23 show the anion-site, p-like 
deep levels in these materials. Here, cation-site means that the 
nearest neighbors of the interstitial are anion atoms, and anion- 
site means that the neighbors of the interstitial are cation atoms. 

Inspecting these figures, we can find some interesting properties 
of the deep levels. First, one can find that cation-site, self- 
interstitial impurities often produce deep levels in the energy band 
gap. In the case of MCT, both self-interstitial impurities Hg and Cd 
produce deep levels well within the band gap. Hg interstitial in 
MZT also produces a deep level within the band gap and Zn 
interstitial in MZT produces a deep level just below the valence 
band edge (this level will move into the band gap after 
considering the effects of lattice relaxation). Hg interstitial and Zn 
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interstitials also produce deep levels in MZS. In contrast with the 
deep levels produced by substitutional impurities, the cation-site 
deep levels produced by interstitials are often acceptor-like. For 
example, in the case of MCT, the deep levels produced by Cu, An, 
Ag, Mg, and Zn interstitials are all acceptor-like. This conclusion 
holds for MZT and MZS. On the other land, the anion-site 
interstitial deep levels can be donor-like or acceptor-like. For 
example, in the case of MCT, N and O interstitials introduce 
acceptor-like deep levels in the band gap, and F and Ar introduce 
donor-like levels in the band gap. 

(b) Effects of lattice relaxation 

The results of the effects of lattice relaxation on interstitial 
deep levels in MCT, MZT, and MZS are summarized in Fig. 
24 to Fig. 29. Most of discussions in section I IF (c) for 
substitutional impurities are also suitable here. For 
example, we can see from these figures that the effects on 
cation-site, s-iike levels are much smaller than those on 

anion-site, p-like levels. An outward relaxation ( a < 0 ) 
moves anion-site, p-like levels toward the conduction band, 

and an inward relaxation (a > 0 ) moves these levels toward 
the valence band. On the other hand, the cation-site, s-like 
deep levels will move the other way under the lattice 
relaxation. 

The effects of lattice relaxation on deep levels produced by 
various interstitial impurities are summarized in Table 2. It 
can be seen that the percentage of the relaxation in all the 
cases is in the range of 0.0% - 16%. N interstitial in MZS 
causes almost no relaxation, and Zn interstitial in MCT will 
introduce only 0.3% relaxation. On the other hand, F 
interstitial in MZT introduces a 16% relaxation. It can be 
also seen that Mg interstitial in MCT, O and Ar interstitial in 
MZT, and Hg and O interstitial in MZS will move out of the 
band gap after considering the effects of lattice relaxation. 
But Zn interstitial in MZT may move into the band gap after 
considering the effects. 

4. Comparison with Experiments 

Many experimental techniques can be used to study the properties 
of defects in narrow gap semiconductors, especially to detect the 
defect energy levels in the band gap. These techniques include 
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photoluminescince (PL), deep level transient spectroscopy (DLTS), 
admittance spectroscopy (AS), thermally stimulated current (TSC), 
magnetic optical spectroscopy (MOS), electron paramagnetic resonance 
(EPR) and so on. Detecting the defect energy levels and identifying them 
are not an easy work to do. For the last two decades, many efforts have 

been made by experimentalists on MCT( c1 ). However, to our knowledge, 
except for one or two cases( c2 ), very little work has been done for 
MZT and MZS. Some experimental data which are firmly identified 
for MCT are shown and compared with our predicted deep levels 
in Table 3. These are the only data available to us at the moment 

It can be seen from Table 3 that our predicted deep levels produced by 
Hg interstitial and Au interstitial are almost the same as the experimental 
data. In other cases, our predicted deep levels are more or less in 
agreement with the data. The differences between the data and the 
predicted levels are less than 0.15 eV. We believe that this discrepancy 
can be explained by our neglecting the effects of charge state splitting in 
the calculations. From this comparison, we have reason to believe that 
our predicted deep levels for MZT and MZS are reliable to an accuracy 
on the order of 0.1 to 0.2 eV. The calculations are a considerable 
improvement over the 0.5 eV uncertainty typical of other results. 

5. Comparison with Other Theories 

(a) Lattice relaxation of substitutional impurities in HgTe 

The problem of lattice relaxation introduced by substitutional 
impurities in MCT has also been studied by other research groups 
using different theories. For example, Berding and Sher use a total 
energy method to calculate the defect formation energy in MCT 
which can be used to predict the relaxation about the defect 
site( c3 L In their approach, defect formation energies are calculated 
from the difference in total energies of the compound with and 
without the defect. They found that for the composition x=0, a Hg 
antisite defect ( in a Te site ) and a Te antisite defect (in a Hg site ) 
in HgTe introduce 0.06 A and 0.2 A outward relaxation, 
respectively. In our calculation we found that the corresponding 
relaxations are 0.05 A and 0.13 A which are reasonably close to 
their results. 

(b) Interstitial impurities in MCT 

The problem of interstitial impurities in MCT is also studied by 
Morgan-Pond, Goettig, and Schick using a total energy 
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method( c4 c6 ). In order to calculate the lattice relaxation induced by 
an interstitials, they use a supercell aproach in which there are 64 
atoms in the bulk crystal supercell and one more for the interstitial. 
They predict that Cd, Hg, and In interstitilas in MCT may have 3% 
to 4% relaxation. From Table 3 we can see that Cd and Hg 
interstitials in MCT are predicted by our calculation to have 4.3% 
and 3.6% respectively, which are in agreement with their 
predictions. Morgan - Pond et al. also calculate deep levels 
produced by Indiun and self interstitials (Hg, Cd, and Te) in MCT 
using a cluster method. They found that a Hg interstitial may form a 
s - like level in the band gap, and an In interstitial may produce a s 
- like level resonant just below the valence band. From Fig. 19 we 
can see that the Hg interstitial is predicted to form a s - like deep 
level at about 0.8 Eg, and a s - like In interstitial deep level is 
predicted to resonate in the valence band. Essentially, our 
predictions agree with theirs. 

ill. Work can be Done Later 

In the present work, we developed a formalism for investigating deep defect states in 
narrow gap semiconductors. With the formalism, we can study the electronic properties 
of either impurities or native defects (including antisite defects, vacancies, and 
interstitials). Specifically, the deep levels produced by these impurities and defects in 
the band gap of these materials can be predicted. The major advantage of this 
formalism is that different types of defects can be easily studied for alloys of various 
compositions of one or more types of materials. In comparison to other calculations, 
the accuracy of our calculations is much improved by including the effects of lattice 
relaxation on deep levels. One can see from Table 3 that the uncertainties of our 
predicted deep levels are less than 0.15 eV in comparison with the experimental data. 

We note that, however, our formalism has neglected the effects of charge state 
splittings on deep levels produced by ionized impurities and defects. This formalism 
can, thus, be used only for neutral impurities and defects. In most cases where data 
are available, the charge state splittings of deep levels are of the order of 0.1 eV to 0.3 
eV. For example, the measured splittings of S, Se, and Te in Si are 0.295, 0.286, and 
0.212 eV( d1 ). Myles predicted that the splittings of deep levels in Hgl-xCdxTe are of 

the order of 0.1 eV to 0.2 eV( d2 ). A sophisticated deep level theory also predicted that 
a neutral Ga vacancy in GaP has its energy level to be 0.2 eV lower than the 

negatively charged vacancy( d3 ). Since the band gap of a narrow gap semiconductor is 
typically 0.1 eV, the effects of charge state splitting are very important in the 
calculations of deep levels. Based on these facts, we have reason to believe that the 
remaining uncertainties in our predicted deep levels can be explained by the effects of 
charge state splitting. The formalism can, thus, be further improved by including these 
effects. 
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We also note that our formalism failed to study the effects of lattice relaxation on deep 
levels produced by vacancies. The reason of this failure is that we use the “ideal 
vacancy model” in the calculations of vacancy deep levels. That is, when we calculate 
vacancy deep levels based on this model, we assume that the diagonal elements of 
the defect potential matrix are infinity. This is equivalent to remove an atom from the 
crystal, but leaving all other atoms in the same position. Thus, the effects of lattice 
relaxation which are represented in the off-diagonal elements of the defect potential 
matrix have no affect on vacancy deep levels. A more sophisticated theory, however, 
indicated that a 10% lattice relaxation may cause a 0.3 eV shift of a deep level with no 
relaxation^ 3 ). It is worthy to note that a Cd vacancy may form a deep level at about 0.4 
eV above the valence band edge in CdTe( d4 ), but our formalism predicted this level at 
about 0.1 eV above the valence band edge in CdTe. We believe that this conflict can 
be explained by the effects of lattice relaxation. Obviously, in order to improve the 
accuracy of the calculations of vacancy deep levels we must construct a better model 
than the “ideal vacancy model”. 

The charge state effects can be included in the calculations of deep levels in an 
approximate manner, using a self - consistent procedure. Following Haldane and 

Anderson( d5 ), the atomic - orbital energies, or ionization potential, in different charge 
states are approximated by 

Esolin }) = E° + ^ n so U ss + ^ n p,, U sp 

c j,a 


Epa(fn}) Ep + ^ ^pa Upp + ^*sa ^sp 

j, a a , (9) 

where a is the spin. Here n S(J and n pa are the occupation numbers ( 0 or 1) of the s 

orbital and p orbital of spin a, respectively. The collective set of occupation number is 
denoted by {n}.The quantities E^ s , E^p, U ss , U sp , and U pp are empirical parameters 

chosen so that they reproduce i) the Hartree - Fock occupied s and p orbital energies 
for the neutral free atom and ii) the experimental s and p ionization potentials. Sankey 

and Dow have determined them for numerous atoms^ 6 ). The diagonal matrix elements 
of the defect potential will now be 


(10a) 


u sa= IW E so( im P) " E sa( host )]’ 
and 

^pa = Ppa[ E pcr(' m P) " E pa(^ os *)]- 
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(10b) 



The spin - orbit occupation number n a ( a = st , si , p| T , p| l , i = x, y, or z) can be 
written as< d7 > 


n a = /_°F(E)D(E)dE, 


(11) 


where F(E) is the Fermi function, D(E) is the local spectral density of states which can 
be written as 


D(E) = <p |p (E) | p> = — (1 /jc) lm<p | G(E) I p> = -(1/«) Im <|i |[ 1 — G0V]-1 1 jj>, (12) 
where G 0 is the host Green’s function and V is the defect potential matrix. 

For a particular impurity, a scheme for obtaining self - consistent solutions is 
implemented as follows. For an input defect potential V, the deep level energy E can 

be calculated by Eq(1). With this E and V, the spin - orbit occupation numbers n are 

oc 

computed from Eq(1 1). The n a are then used to compute the spin - orbit energies Eso 

and Epofrom Eq(9). These energies are then used to compute new potential matrix 
elements from Eq(10 ). Finally, the resulting V is used in Eq(1) to calculate a new deep 
level energy E. This procedure is repeated iteratively until self - consistency is 
obtained. 


The effects of lattice relaxation on vacancy deep levels can be considered following 
Singh and Madhukar’s model( d8 ). To simplify the problem, let’s consider a cluster with 
two atoms per unit cell and one s orbital per atom. We use the subscripts 0 and 1 to 
denote these atoms. The Hamiltonian of the perfect crystal can then be written as 


H, 


■-S' 




+ Si>V<j| 
>.j 


(13) 


where p takes 0 and 1 , and i and j are nearest neighbors. If a defect is at 0 site, then 
the interactions describing the defect are given by 


< 0, s 0 I H 1 0, s 0 > = Eg 

<0, s 0 I H 1 1, Si > = V 01 

< 1, s-j I H 1 1, s-j > = e^ 

< 1, s-j I H 1 2, Sq > = V-j2» (14) 


where H = H 0 + V is the Hamiltonian including the defect, and 2 denotes the 
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neighboring unit cell. Following Singh and Madhukar, the diagonal element of the 
Green's function matrix which includes a defect is given by 


Q oo = 


(E-s 0 ) 


4V 


0 1 


3V. 




0 VT 0 


where T 0 is given by 


(15) 


T 0 ~v 


(E ~ e 0 ) - 


3/ 2 

1 

c „ 4V 2 

1 

E -£ - 

1 

0 

G oo _ 


(16) 


For an ideal defect at site 0, we have V 12 = V, and £i = Therefore, Eq(16) becomes 


G oo = 



(17) 


3n the “ideal vacancy” model, e— > oo , G 00 , thus, is zero. The deep level vacancy can 

then be determined. But using this model, the vacancy deep levels are determined by 
the relation 


( E ‘ e i) " [(e - £i) + a] — = 0 , 

J v 

(18) 

where A = 4V 2 0 f G 00 - The effects of lattice relaxation are included in V 0 -j and V^. so 
that Eq(18) gives these effects on vacancy deep levels. 
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Tablet: Effects of lattice relaxation on substitutional deep levels. 


System 

d H (A) 

d,(A) 

Ax=d| - d H 

% of relaxation 

deep levels 3 
(no relax) (relaxed) 

MCT: Cdy e b 

2.8 

2.97 

0.17 

6.1 

C.B.° 

C.B. 

MCT: Mgy e 

2.8 

2.96 

0.16 

5.7 

°- 3E g d 

a5E g 

MCT: Znj e 

2.8 

2.98 

0.18 

6.4 

°.5E g 

0. 9E g 

MCT: Hgy e b 

2.8 

2.97 

0.17 

6.1 

C.B. 

C.B. 

MZT: 

2.74 

2.61 

-0.13 

4.7 

0.98Eg 

0.7Eg 

O.SEg 

MZT: Mgj- e 

2.74 

2.60 

-0.14 

5.1 

0.5E a 

MZT: C dj- e 

2.74 

2.60 

-0.14 

5.1 

C.B. 

C.B. 

MZT: Hgy e b 

2.74 

2.60 

-0.14 

5.1 

C.B. 

C.B. 

MZS: Beg e 

2.62 

2.74 

0.12 

4.6 

O.SEg 

0.6Ea 

MZS: TI Se 

2.62 

2.77 

0.15 

5.7 


0AE° 

0.SE 

MZS: lng e 

2.62 

2.77 

0.15 

5.7 

-'•<4 

MZS: Hg Se b 

2.62 

2.73 

0.11 

4.2 

C.B. 

C7 

C.B. 

MZS: Zn Se b 

2.62 

2.73 

0.11 

4.2 

C.B. 

C.B. 


a. AH deep levels in this table are anion site, p-like deep levels. 

b. All these are antisite impurities. 

c. C.B. means conduction band. 

d. Eg means energy band gap. 
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Table2: Effects of lattice relaxation on interstitial deep levels 


System 

d H (A) 

d|(A) 

Ax=dj - d H 

% of relaxation 

deep levels 3 
(no relax) (relaxed) 

MCT: Cc*> 

2.8 

2.94 

0.12 

4.3 

0.74E g d 

a77E g 

C.B. C 

MCT: Mg 

2.8 

3.10 

0.24 

8.6 

0.98E a 

0 - 1E g 9 

MCT: Zn 

2.8 

2.81 

0.01 

0.3 

0,1 E g 

MCT: Hep 

2.8 

2.90 

0.10 

3.6 

0.82 E g 

0.83E g 
0.91 Eg 

MCT: Au 

2.8 

2.90 

0.10 

3.6 

0 . 9Eg 

MCT: Ag 

2.8 

2.91 

0.11 

3.9 

0 . 87Eg 

0.88E g 

MCT: Cu 

2.8 

2.91 

0.11 

3.9 

°- 2E g 

°* 22E g 

MCT: T& 

2.8 

2.95 

0.15 

5.4 

C.B. 

C.B. 

MCT: N 

2.8 

2.56 

-0.24 

8.6 

0 . 9Eg 

°- 7E g 

MCT: 0 

2.8 

2.54 

-0.26 

9.3 

0 . 8Eg 

°- 6E g 

MCT: Ar 

2.8 

2.59 

- 0.21 

7.5 

0.82E g 

0.62E g 

MCT: F 

2.8 

2.52 

-0.28 

10 

0,44E g 

a24E g 

MZT: ZrP 

2.74 

2.92 

0.18 

6.6 

V.B. C 

0.04E g 

MZT: Mg 

2.74 

2.80 

0.06 

2.2 

O.OEg 

0.6Eg 

0.91 Eg 

MZT: Cd 

2.74 

2.82 

0.08 

2.9 

0.61 Eg 

MZT: Hep 

2.74 

2.81 

0.07 

2.5 

0.68E g 

0 . 69Eg 

MZT: Au 

2.74 

2.80 

0.06 

2.2 

°- 74E g 

0.55E g 

MZT: Ag 

2.74 

2.80 

0.06 

2.2 

0 . 70Eg 

°- 71E g 

MZT: Cu 

2.74 

2.88 

0.14 

5.1 

0.1 8E g 

0.30E g 

MZT: Td> 

2.74 

2.35 

-0.41 

15 

C.B. 

C.B. 

MZT: N 

2.74 

3.10 

0.36 

13 

C.B. 

C.B. 

MZT: 0 

2.74 

3.14 

0.40 

14.6 

0 . 8Eg 

°- 9E g 

C.B. 

MZT: Ar 

2.74 

3.15 

0.41 

15 

C.B. 

MZT: F 

2.74 

3.18 

0.44 

16 

a3E g 

O.SBg 
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Table 2 : (continued) 


MZS: Hg b 

2.62 

2.71 

0.09 

3.4 

0.98E g 

C.B. 

MZS: Zn b 

2.62 

2.78 

0.16 

6.1 

0.04Eg 

°- 14E g 

MZS: Cd 

2.62 

2.72 

0.10 

3.8 

0.86E g 

0.88E g 

MZS: Cu 

2.62 

2.75 

0.13 

5.0 

°.24 E 

0.34E a 

y 

MZS: Ag 

2.62 

2.71 

0.09 

3.4 

C.B. 

C.B. 

MZS: Se b 

2.62 

2.58 

-0.04 

1.5 

C.B. 

C.B. 

MZS: N 

2.62 

2.61 

-0.01 

0.0 

C.B. 

C.B. 

MZS: O 

2.62 

2.37 

-0.25 

9.5 

0.98E g 

C.B. 

MZS: F 

2.62 

2.51 

-0.11 

4.2 

0.50E g 

0.40E, 


a. All deep levels in this table are either anion site, p-like deep levels or cation site, 
s-like deep levels. All interstitials are on the preferred positions. 

b. All these are self-interstitials. 

c. C.B. means conduction band edge, and V.B. means valence band edge. 

d. Eg means energy band gap. 
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Table 3: Comparison of the predicted deep levels from the present theory 
with experimental data 


system 

deep levels 


experiment 3 

theory 

Hgj b (donor) 

O.SEgC 

0.83E g 

Tej b (recom) d 

°.5Eg 

E c e + o. 1 eV f 

Te Hg ( recom ) 

O.SEg 

E c + 0.06eV 

Cu Te (donor) 

Q.8Eg 

E c + 0.06eV 

Au-j- e (donor) 

0.8Eg 

E c + 0.04eV 

Auj b (donor) 

°.8Eg 

0.91 Eg 


a. Data obtained from References el - e3. 

b. These are interstitial impurities. 

c. Eg means energy band gap. 

d. Recombination center. 

e. E c means conduction band edge. 

f. The difference between data and theory result can be explained by the effects 
of charged-state splitting. 

g. The x values for the experimental data here are in the range 0.2 < x < 0.3. for 
our calculations, x= 0.22. As seen by Figs. 1 and 2, the change in the energy 
levels for 0.2 < x < 0.3 is very small and can be ignored for rough comparisons 
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JV. 


Results for Shallow Levels in II-VI Materials 


In the mid year report* 31 ), we already gave a fairly complete discussion of 
the basic theory of shallow impurity levels in II-VI materials. In the very 
lowest approximation we obtain the simple hydrogenic model (with an 
appropriately scaled Bohr radius and Rydberg energy) but of course this 
is not a good approximation for our materials. To get better results, we 
must include the effects of the following. 

1. Non-parabolic Bands 

2. Screening 

3. Site Dependence 

4. Central Cell Corrections 

5. Alloy Disorder 

6. Lattice Coupling 

7. Lattice Relaxation 

8. Impurity Interactions 

9. Band Coupling and Degeneracy 

From previous work to first order* 31 ), the Hamiltonian that corrects for 

non-parabolic behavior is 

H 1 = ap 2 + bp z <|> p z + cp + <J>p_ + dp.<|>p + n 


where P is the momentum operator and <|> is the hydrogenic potential. 
The coefficients a, b, c, and d are given in the mid year report. 

Using appropriate values for the parameters for MCT and MZT we can 
evaluate E gap , E donor , ground state and the ground state donor wave 

function \\ r 0 . To first order 


AE = \\|f 


H 




(20) 


Our results are given in Table IV. The results predict an approximate 
10% downward shift in binding energy due to the effect of non- 
parabolicity. Of course, many other factors come into play as noted 
above. Perhaps the most interesting part of this is that we predict an 
increase in binding due to non-parabolicity. 
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TABLE 4 



MCT 

MZT 


1.290 x10 19 eVcm 

1.238X10 19 eVcm 

p 

erg sec 

erg sec 

X 

.19 

.13 

K 

17.56 

18.35 

A 

1.05 eV 

.9 eV 

T 

300 K 

300 K 

Eg 

.16 eV 

0.141 eV 

Donor 



“-Ground 

- 0.014 eV 

- 0.01443 eV 

AE 

-0.00146 eV 

-0.00211 eV 

AE 



EGap 

.9% 

1.5% 

AE 



— Donor 



Ground 

- 10.4% 

- 14.6% 

T (same p, x, k, A) 

77K 

77K 

Eg 

.0841 eV 

0.0800 eV 

_ Donor 



•-Ground 

- 0.008842 eV 

- 0.00828 eV 

AE 

- .001268 eV 

-0.001 16 eV 

AE 



EGap 

1.5% 

1.45% 


Results for shallow levels in ll-VI Materials 
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MOT 


MZT 


AE 


E 


Donor 

Ground 


1 .45% 


- 14 % 


E g (same p, x, k, A) 

— Donor 
^Ground 


AE 




Donor 

Ground 


.1 eV 

- 0.010082 eV 

- 0.001352 eV 


1.4% 


- 13.4% 


.1 eV 


- 0.009785 eV 
-0.0012464 eV 


1 .25% 
- 12.7% 
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(antisite) 


Fig. 2 Substitutional Deep Levels in MCT as 
a Function of x (anion - site, p - like) 





band gap energy (eV) 


0.4 



Fig. 5 Substitutional Deep Levels in MZS as 

a Function of x (anion - site, p - like) 


(antisite) 
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band energy (eV) 



Fig. 8 Substitutional Deep Levels in MZT for 
Eg = 0.1 eV (cation - site, s - like) 
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band energy (eV) 



Fig. 9 Substitutional Di 
Eg = 0.1 eV 0 
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band energy (eV) 



impurity potential (eV) 


Fig. 11 Substitutional Deep Levels in MCT for 
Eg = 0.1 eV with Lattice Relaxation 
(cation-site, s - like) 
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band energy (eV) 



impurity potential (eV) 


Fig.12 Substitutionai Deep Levels in MCT 

for Eg = 0.1 eV with Lattice Relaxation 
(anion - site, p - like) 
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band energy (eV) 



impurity potential (eV) 


Fig. 15 Substitutional Deep Levels in MZS for 
Eg = 0.1 eV with Lattice Relaxation 
(anion-site, p - like) 
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band energy (eV) 



Fig. 17 Vacancy Deep Levels (anion - site, 

p - like) 
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band energy (eV) 



Fig. 18 Interstitial Deep Levels in MCT for 
Eg = 0.1 eV (cation - site, s - like) 
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band energy (eV) 


0.2 


1 


0.1 H 


C.B. 


N 


A r 


0.0 H 


V.B. 


-o.i -4 


J 


Fig. 19 interstitial Deep Levels in MCT for 

Eg = 0.1 eV (anion - site, p - like) 
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band energy (eV) 





band energy (eV) 



Fig. 21 Interstitial Deep Levels in MZT for 

Eg = 0.1 eV (anion - site, p - like) 
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band energy (eV) 



Fir. 22 Interstitial Deep Levels in MZS for 
Eg = 0.1 eV (cation - site, s - like) 



band energy (eV) 


0.2 



Fig. 23 Interstitial Deep Levels in MZS for 

Eg = 0.1 eV (anion - site, p - like) 
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band energy (eV) 



impurity potentia! (eV) 


Fig. 24 Interstitial Deep Levels in MCT for 

Eg = 0.1 eV with Lattice Relaxation 
(cation - site, s - like) 
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band energy (eV) 


impurity energy (eV) 



Fig. 25 


Interstitial Deep Levels in MCT with 
Lattice relaxation (anion - site, 
p - like) 


band energy (eV) 



Fig. 26 interstitial Deep Levels in MZT for 
Eg = 0.1 eV with Lattice Relaxation 
(cation - site, s - like) 


48 


impurity potential (eV) 


interstitial Deep Levels in MZT for 
Eg = 0.1 eV with Lattice Relaxation 
(anion - site, p - like) 



band energy (eV) 



impurity potential (eV) 


Fig. 28 interstitial Deep Levels in MZS for 
Eg = 0.1 eV with Lattice Relaxation 
(cation - site, s - like) 
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band energy (eV) 



impurity potential (eV) 


Fig. 29 interstitial Deep Levels in MZS for 

Eg = 0.1 eV with Lattice Relaxation 
(anion - site, p - like) 
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Space Grown Crystals 


Jim Patterson, Sandi Billings, and 
JimMantovani 

Tie exploitation of space began with the 
inch of Sputnik in 1957, which started ' 
; race to the moon and the race for 
liars. In the name of national pride, tax 
pport was rallied not just for exploring 
e unknown, but also because of the 
omise that scientific returns will make 
oney through the development of new 
oducts and, in the process, will raise the 
andard of living. Of course, space has 
;en used for military purposes and for 
irious benign remote sensing purposes, 
at the idea of manufacturing something 
i space is particularly attractive, promis- 
lg a return on the dollar. After more than 
5 years of promises, we may have found 
space-based industry that will deliver. 
The growing of crystals in space meets 
aany of the criteria for space manufac- 


CRYSTAL 

CdZnTe 

HgCdTe 
HgZnTe 
GaAs (doped) 

InSb 

PbSnTe 

Protein 


The basic reason 
defect free, homogeneous crystals. 


ture. The product — nearly defect-free 
crystals — is extremely valuable, com- 
pact, and light, and made much better 
than on Earth. The growth of certain 
kinds of crystals seems to benefit from the 
microgravity environment of space^. 2) 
The profitability of manufacturing crys- 
tals in space, however, remains an open 
question. 

In this article we first review what crys- 
tals are and how they are grown, and then 
we discuss the concept of microgravity 
and why it might be an attractive environ- 
ment for growing crystals. Before discuss- 
ing the Space Station and the Space 
Shuttle, we will review some earlier tech- 
niques for growing crystals in 
microgravity. Finally, we will discuss 
some recent and proposed techniques for 
growing crystals in space. 

Most people think they know what a 
crystal is — after all, it’s as clear as crys- 


USE 

Material Substrate 
for HgCdTe 

Infrared detection 

Infrared detection 

High Speed IC’s and 
Lasers 

Various electronic devices 

Broadband and 
Infrared detection 


tal, isn’t it? Actually, crystal glass isn’t 
crystalline at all; it is a special kind of 
amorphous or random structure. A crys- 
talline structure has a regular array of 
atoms, as regular as the arrangement- of 
black and white squares on a checker- 
board. This perhaps can give us a clue as 
to why it is difficult to grow perfect crys- 
tals on Earth under the influence of grav- 
ity. 

Crystal Defects 

No crystal is perfect. Some defects may 
be desirable, such as the donors and ac- 
ceptors introduced into semiconductors to 
change their electrical conductivity. Even 
when desirable, the number and types of 


About the authors: 
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for growing all crystals in space is to produce large; 
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defects must be controlled, or the crystal 
won’t behave as desired. The kinds of 
defects that occur are point, line, and 
plane defects. 

Point defects come in many varieties: 
vacancies, self and foreign interstitials, 
substitutional impurities, and even 
antisite defects. Several point defects are 
illustrated in Figure 1. Antisite defects 
only arise in compound materials with at 
least two kinds of atoms. A typical com- 



Figure 1. Point defects: (a) interstitial 
atom— self or foreign, (b) substitutional im- 
purity, (c) vacancy-absence of atom. 


pound ciystal of interest is shown in Fig- 
ure 2. An antisite defect occurs if an A 
atom is found on a B site or vice versa. 
Several of the crystals we will discuss are 
compound. 

Line defects come in two types: edge 
and screw dislocations, which are shown 
in Figure 3. An edge dislocation is simply 
a missing half-plane of atoms. A screw 
dislocation, occurring when one part of a 
crystal is pushed up relative to other parts, 
can be visualized as cutting the pages of a 
book halfway across at the middle. Then 
the lower part of page 1 can be “taped” to 
the upper part of page 2 and similarly for 
the following pages. This would form a 
“spiral staircase.” Except for the “tap- 
ing,” the idea is shown in Figure 3b. 

Most crystals are really just collections 
of many small, individual crystals called 
grains. Grains come together at grain 
boundaries, which are examples of plane 
defects. For many purposes it is desirable 
not to have grains, but rather a large 
single crystal. 

Crystal Growth 0,4 * 

One of the ways of growing crystals on 
Earth is the Bridgman technique, which 



Figure 2. ZnS structure. A atoms (e.g., Zn) 
are located at sites to (0, 0, 0). B atoms (e.g., 
S) are located at sites equivalent to (a/4, a/4, 
a/4), where a is the length of the side of the 
unit cell shown. 

also works in microgravity. In a common 
Bridgman process, shown in Figure 4, the 
material is melted in a stationary vertical 
column which tapers to a point at the 
bottom where there is an unmelted 
“seed” crystal. As the furnace moves up, 
the crystal resolidifies from the seed. The 
steps in the crystallization process are dif- 
fusion of molecules to the surface of the 
crystal (the solid melt interface); diffusion 
of molecules over the surface; binding of 
the molecules at the surface; and diffusion 
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leat away from the surface of crystalli- 
on. When not too far from equilibrium 
ditions, screw dislocations can be an 
ortant aid to growth, because the at- 
; have a higher probability of sticking 
ie step caused by the screw. Very large 
artures from equilibrium can cause 
ilar growth to be unstable; it is then 
aced by dendritic growth, forming 

thin crystals. 

addition to the Bridgman technique, 
e commonly used melt growth tech- 
ies are floating zone, growth from so- 
>n, and growth from vapor. The float- 
1 method, which has some similarities 
ie Bridgman technique, uses a rod of 
material (not single crystalline) and a 
le crystal rod of the same material, 
.rods are joined end to end and heated 
1 a melt is formed. This melt is moved 
the raw materials, and the crystal is 
yn from the seed. In microgravity 
' is no distortion of the molten zone, 
convection can still be a problem. In 
t crystal growth the sample is heated 
vaporized. It then flows to the cooler 
>n and is deposited on the substrate, 
vth from solution is important prima- 
n protein crystal growth, 
ravity causes convection, buoyancy, 
hydrostatic pressure. Convection, 
cularly turbulent convection, can 
' many kinds of defects to be intro- 
i into the crystal. Without convec- 
crystal growth can occur under diflfu- 
controlled conditions, leading to im- 
cd chemical homogeneity and struc- 
perfection. It is also possible to grow 
i larger single crystals than could be 
n under the stress caused by gravity, 
desired goal is to grow large single 
ils which are nearly perfect. Once 



e 4. Crystal growth by the Bridgman 
ique. T a is the melting temperature of 
ystaL 
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Figure 3. Dislocations: (a) edge and (b) screw 



this has been accomplished, controlled 
amounts of defects can be introduced as 
desired. In the absence of convection and 
other stresses introduced by gravity, 
nearly perfect single crystals are much 
more likely to form. 

In microgravity, some effects that are 
totally masked by gravity-driven causes 
show up nicely, such as Marangoni con- 
vection. In a melt, the surface tension may 
vary with temperature, and in the absence 
of gravity this can be enough to cause 
convection— Marangoni convection, 
which may cause striations in crystals 
grown in space. 

Microgravity 

First we should talk about what 
microgravity is and what it is not. It is not 
the absence of gravity, or even a region 
where gravity is very small. Unless one 
goes very far from massive bodies, this is 
impossible. Even at a Shuttle orbit of 300 
km above the Earth, the force of gravity is 
91% the value experienced on the Earth. 
Microgravity is accomplished by having 
almost all the effects of gravity canceled 
out by free fall. 

Newton himself understood the prin- 
ciple, illustrated in Figure 5. If one 
mounts a cannon on a large mountain on 
an otherwise flat Earth and fires the can- 
non ball horizontally, it will land some 
distance from the base of the mountain. 
Adding more powder will cause the ball to 
go further. Finally a point will be reached 
when the ball falls exactly the same 
amount that the Earth curves. The ball 
will then be in free fall, and in orbit. The 
effects of gravity for objects inside the ball 
will be very small. In an orbiting satellite 
there will be exactly one surface where the 


effects of gravity are negligible and at other 
places inside one has “microgravity.” 
There are many ways to produce 
microgravity; all you have to do is arrange 
to be in free fail. (2) Drop towers and drop 
tubes offer two ways of accomplishing 
this. The first commercial use of micro- 
gravity was probably the drop tower used 
in 1785 in England to make spherical lead 
balls. Marshall Space Flight Center has 
both a drop tower and a drop tube 100 
meters high - this alone allows one to have 
free fall, or microgravity, for about 5 sec- 
onds. In a drop tower the entire experi- 
mental package is dropped. For crystal 
growth experiments, this means the fur- 
nace as well as the instrumentation and 
the specimen are all placed in a special 
canister and dropped. In a drop tube, there 
is an enclosure in which, for example, 
only the molten sample would be dropped. 
Special aerodynamic design, vacuum, or 
other means is used to reduce air drag and 
hence obtain real free fall. For slightly 
longer times (20 seconds or so), the KC 
135 aircraft can be put into a parabolic 
path to produce microgravity. Extending 
this idea, rockets have been used to pro- 
duce microgravity for periods of about 
400 seconds. 

For long-duration experiments, the ex- 
perimental package needs to be in orbit. 
Such research began with SKYLAB in 
1973 and has continued with the Space 
Shuttle and the Space Lab in 1983. We 
now look forward to Space Station Free- 
dom and to spending very long times in 
space. This is important for crystal 
growth, because some crystals grow so 
slowly that the process is about as exciting 
as grass growing. And the longer that we 
have to grow crystals in space, the larger 
we can grow them. 
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Crystals and Crystal Growth Furnaces 
for Space 

Many crystals are being considered for 
growth in space. Among the more inter- 
esting of them are the narrow gap semi- 
conductor alloys. Mercury Cadmium Tel- 
luride (MCT), Mercury Zinc Telluride 
(MZT), and Mercury Zinc Selenide 
(MZS) (5) . At Florida Tech in Melbourne, 
FL, authors Jim Patterson and Jim 
Mantovani are studying these materials 
under NASA research contracts. One use 
of these semiconductor crystals of interest 
to NASA is making infrared detec - 
tors (8 ’ 7,8) . 

Because of its stability with respect to 
mercury, MZT appears to be a strong 
competitor to MCT. It is very hard to 
produce homogeneous bulk crystals of 
MZT on Earth due to convection and 
composition segregation driven by grav- 
ity. We also need to point out that these 
compound systems have phase diagrams 
which must be considered in predicting 
the nature of the crystal formed. (5) For 
MCT we see in Figure 6 that the fraction 
of CdTe which solidifies depends on the 
temperature. 

Scientists and engineers are presently 


hard at work producing crystal growth 
furnaces for use in microgravity. A typical 
one, used by Dr. S.L. Lehoczky of 
Marshall Space Flight Center, (3) is shown 
in Figure 7. It has three major elements: 
(1) an avionics subsystem (the electronics 
to control its operation); (2) the environ- 
mental control system, obviously used to 
produce the proper environment in either 
fluid or vacuum; and (3) the furnace, 
which consists of several parts. Typically 
the furnace is set up to have several heat- 
ing zones which produce a temperature 
gradient. It is found that it is better to 
move the furnace itself rather than the 
crystal, so as not to introduce instabilities 
into the growing process. Even so, move- 
ments in the Shuttle can produce g-jitter. 

Besides MCT, MZT, and MZS, there 
are several other crystals being considered 
for growth in microgravity, as listed in 
Table I. These include InSb, GaAs, and 
protein crystals. However, the whole sub- 
ject of crystal -growing is far more compli- 
cated than many would imagine. (9 I0) De- 
spite much initial excitement, the viability 
for commercial growth of crystals in space 
is an open question, which a number of 
Brevard industries and universities will 
continue to explore. 
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NARROW GAP SEMICONDUCTORS 


by 

J.D. Patterson 


Florida Institute of Technology 
Melbourne, FL 32901-6988 
U.S.A. 


Introduction 

Not all semiconductors are alike. Narrow gap semiconductors (NGS)* 1 * 7 ) are often 
defined as having an energy gap, as measured from the top of the valence band to the 
bottom of the conduction band, of less than or about 0.5 eV. This is to be contrasted 
with the energy gap in Si at room temperature of about 1.1 eV. NGS can be regarded 
as being intermediate between covalently bonded solids and metals. Because of the 
narrow gap between them, there is a strong interaction between the conduction and 
valence bands and the energy dispersion relation becomes non-parabolic (E not 
proportional to k 2 ). Also, since the effective mass of carriers depends on their kinetic 
energy at the Fermi surface, the mobility may depend on the carrier concentration. 

NGS are sensitive to external influences such as temperature, external magnetic 
fields, and strain, features which make them strong candidates for sensors of many 
kinds. Besides their principal use as infrared detectors and emitters, NGS have 
application or potential application in magnetoresistive, Hall effect, and thermoelectric 
devices, and in semiconductor strain gauges. IR detectors, of course, have a large 
variety of uses including military (night vision and detecting the plume of ballistic 
missiles) and peaceful (remote sensing and security systems). 

As has been observed* 8 ), narrow gap semiconductors have been in use for a long 
time. Their first recorded use was as eye liners for the ladies of early Egypt. For this, 
crushed Galena (PbS, lead sulphide) was used. A still later use of Galena was as a 
crystal in early crystal detector radios. PbS is a member of the IV-VI lead salt family 
which also includes PbSe and PbTe. They have been used as infrared (IR) detectors, 
but their main use nowadays is for lasers. By using several compounds (e.g. PbSnTe, 
PbSe and PbEuSeTe), the wavelengths from 3 to 30 nanometers can be covered. 

More recently discussed materials, such as the lll-V family of semiconductors, have 
found a good deal of application (see Table I). A particularly important member of this 
group is InSb a NGS first discussed in detail by Welker* 9 ) in 1952. Later E.O. Kane* 10 ) 
gave a fairly simple discussion of its band structure. InSb is still an important material 
for use in IR detectors. 
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The materials that we are particularly concerned with in this review are the ll-VI 
compound semiconductors which include mercury cadmium telluride (MCT) and 
mercury zinc telluride (MZT) (see Table II). Lawson* 11 ) was the first to do detailed work 
with the MCT alloys. He showed how one could continuously vary the energy gap by 
varying the ratio of the concentration of Cd to Hg. He thus developed materials which 
operate as IR detectors in the 10 micrometer range of wavelength and at 77K. Several 
other NGS (e.g. Bi 2 T%, Te, and Mg^Sn) are mentioned in the paper by R.S. Ailgaier in ' 
the collection edited by Seiler and Littler* 2 ). 

General Remarks 

In this paper we will focus our attention on the ll-VI compound narrow gap 
semiconductors* 12 ' 14 ). Figs. 1-3 show the zinc blende crystal structure of the materials, 
their first Brillouin zone, and an example of their calculated band structure. ZnTe and 
CdTe are both direct band gap semiconductors. Fig. (4a) shows the energy gap of 
MCT vs. x, the fraction of Cd, at room temperature (upper curve) and at 77K (lower 
curve). One reason MCT is so important is that its energy gap E g can be tuned by 
varying x as we have just shown. In practice, the calculation of the dependence of E g 
on composition and temperature is very difficult and so one often resorts to obtaining 
Eg (x,T) through curve fitting. A good empirical approximation for the energy gap in 
Hg.|_ x Cd x Te is* 15 ) 

Eg (x,T) = -0.302 + 1.93x + (5.35)( 10' 4 )(1-2x)T - 0.81 x 2 + 0.832x 3 (eV) 

for 0 <x < 1 and 4.2 <T< 300K. (1) 

Fig. (4b) shows the energy gap for several materials (at room temperature) and the 
corresponding inverse cutoff wavelength (X c ) given by 



for X c in micrometers and E g in eV. Photons of energy less than Eg are not absorbed 
and X c \s the corresponding wavelength above which no absorption occurs. 

Fig. 5 sketches the energy vs. k dependence for the conduction band, heavy and light 
hole bands, and the split off bands as determined by the model of Kane* 10 ) originally 
for InSb. For further details see the appendix. 

The ll-VI semiconductors, that we are particularly concerned with here, are important 
for several reasons, as indicated in a summary of the properties of the most important 
member of this group (MCT) given in Table 3. Some useful properties of elements 

2 



involved in compound MCT semiconductors are given in Table 4. We go on to discuss 
two particularly important compound semiconductors MCT and MZT. 

MCT and MZT 

We will discuss in some detail two properties which are subjects of current research, 
namely electron mobility and deep defects. * - 

Mobility 

High mobility materials make good IR detectors, among other things, hence the theory 
of mobility is of practical importance. The theory for compound ll-Vi semiconductors is 
more complicated than the theory of simple semiconductors. These complications 
arise due to the non-parabolic nature of the energy band and the highly ionic 
character of the material which makes polar interactions of the electrons with the 
longitudinal phonons play an important role at “high” temperature (high being room 
temperature for many of these materials). In this interaction, the phonon energy is 
comparable to the electron energy. The scattering is inelastic, the relaxation time 
approximation is inappropriate, and a variational principle has to be used to solve the 
non-iterated Boltzmann equation. Another complication is due to narrowness of the 
band gap which implies a non-negligible number of thermally excited electrons, so 
screening has to be taken into consideration. 

The collisions considered in our calculation involve: ionized impurities, holes, 
compositional disorder, acoustic phonons and optical phonons. At low (typically about 
50K) temperature for n-type materials, the scattering is mainly due to ionized 
impurities. This scattering process decreases with increasing temperature, since the 
charged impurity atoms do not scatter high-energy electrons very effectively. The 
theory for hole scattering is somewhat similar to that for ionized impurities. At high 
(typically room) temperature, the scattering is mainly due to longitudinal optic 
phonons. 

Since the two atoms in the unit cell are different, the longitudinal optic mode of the' 
lattice vibrations polarizes the crystal and an electric field is produced. The Frohlich 
perturbation potential is used to describe the interaction of the electrons with the 
polarization field. For the acoustic phonons one uses the deformation potential which 
arises from changes in the energy gap as a result of changes in the spacing of the 
atoms. In a compound semiconductor, like HgCdTe, electrons suffer an additional 
scattering due to the random locations of Hg and Cd atoms. This is termed 
compositional disorder scattering. With known band structure and interactions, the 
Boltzmann transport equation can be solved for the electron energy distribution 
function /, from which the mobility can be obtained. If / is known, we can calculate the 
conductivity by first calculating the electron current density J x , 
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Jx= eV x /D(E)dE, 

J (3) 

where D (E) is the density of states per unit volume and V x is the velocity of the charge 

carriers in the direction of the impressed electric field and current flow. Hence the 
electron mobility can be calculated from 



where n is the number of electrons per unit volume and e is their charge. 
Typical results, obtained by this procedure, are shown in Fig. 6( 16 ). 



(4) 


Defects control the properties of semiconductors. The effects of both shallow and 
deep defects (see Table 5) are considered here. Present understanding of shallow 
defects (which supply carriers) in semiconductors originated with the classic work of 
Kohn and LuttingeH 17 ). The study of deep defects has evolved more slowly but goes 
back to the fundamental paper of Slater and Kosterf 17 ). Deep defects are important for 
both fundamental and industrial research because of their important role in 
recombination of charge carriers (electrons and holes). The study of defects in narrow 
gap semiconductors represents a much less clear story than in semiconductors like 
Si and GaAs. Here, both shallow defects (caused by long range potentials) and deep 
defects (from short range potentials) are far from being completely understood. All 
defects we study are point defects. Until one can identify and characterize all 
appropriate defects, it is impossible to make clear statements about the quality of a 
crystal. A full characterization of the crystal perfection involves not only a theoretical 
study of all common defects but also experimental verification of models for them. The 
simplest model of a shallow defect in NGS assumes the defect is fully ionized but there 
is not a lot one can do with such a model since most of the physics has been 
approximated away. Deep defects may produce energy levels near the center or near 
the band edge< 17 ). Neglect of the long range component (Coulomb tail) for deep 
defects near the band edge is not a good approximation, since it is precisely here that 
small Coulomb tail effects can be important. The Green’s Function model, however, 
gives an accurate estimate of energy trends for deep defects with energies well in the 
gap. The model uses the customary Green’s function^ 7 ), 

G* O) (E) = (E-H 0 | 1 (5) . 

where H 0 is the Hamiltonian of the perfect crystal, and one can show directly from 
Schrodinger s equations that the defect energies E are determined by 
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det (l - G°(E) v)= 0 


( 6 ) 

where V is the defect potential. The Green’s Function can be constructed once band 
structure calculations are available. With a judicious choice of V, this technique has 
been found to produce very good results(^), particularly as regards to chemical trends.' 

There are several kinds of defects as shown in Fig. 7 and Table 6. We study 
substitutional (including antisite, e.g. a Te on a Hg site in MCT), interstitial and ideal 
vacancy defects (see Fig. 7 and Table 7). Vacancies and self interstitials are the most 
common native defects. The ideal vacancy model assumes that a vacancy is formed 
by removing an atom from the crystal, but leaving all other atoms in the same position. 
This is not a very good model and results for this case will not be presented here. For 
substitutional and interstitial impurities, the effects of relaxation of the atoms 
neighboring the substitutional or interstitial impurity are included. Future 
improvements should include charge state interactions (where the impurity is not 
neutral) and dealing with non ideal vacancies in which relaxation is included. 

In considering interstitials, one should note that tetrahedral symmetry for interstitials is 
different from tetrahedral symmetry for substitutional defects because nearest 
neighbors are closer for interstitials than for substitutional impurities. For deep defects 
some results, obtained by us, for the substitutional and interstitial cases are shown in 
Fig. 8a-8d. In all cases we have chosen x so the bandgap is exactly 0.1 eV. For MCT 

this works out to be x = 0.22 at T = OK and for MZT this is x = 0. 1 5 at T = OK. Cation 
sites are mercury, cadmium or zinc sites. Anion sites are tellurium sites. By group 
theory the important levels in tetrahedral symmetry are predominantly s or p like 
(technically A ^ or T 2 ). For the figures, the calculation has been done for each of the 

atoms indicated and after relaxation of neighboring atomic positions, the defect energy 
level comes out as indicated. 

The Future 

An excellent discussion of future applications has appeared. See “New developments 
and trends in the physics and technology of narrow-gap semiconductors,” by R.A. 
Stradling in the book edited by Seiler and Littler**). The future for NGS seems to be 
concerned with applications using new materials and structures* 18 ): multilayers, 
superlattices, multiquantum wells, resonant tunneling devices, strained layers and 
magnetic semiconductors. The surface has barely been scratched as to utilization of 
these new ideas. One of the areas of greatest activity at present is in making IR 
devices from superlattice-like structures. One advantage of superiattices is that one 
can separate the ionized donors (which also scatter carriers) from the region where 
the carriers move. Superiattices also produce “mini bands” which allow us to taylor 
the effective energy gap. Another recent idea is strained layer superiattices. Here a 
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small lattice mismatch between layers can be used to create a strain which modifies 
the electronic properties in a desirable way. Another area of interest is in using 
superlattice type structures to create solar cells of improved efficiencies. Perhaps the 
easiest way to understand how this could occur is to consider layers of decreasing 
band gap. The first layer could remove lower energies in the solar spectrum, the next 
layer could be designed to remove higher energies and so forth, until the whole 
spectrum is efficiently utilized. Experimental techniques of great present and future ' 
importance in order of decreasing expense include MBE (molecular beam epitaxy), 
MOCVD (metalo-organic chemical vapor deposition), and electrodeposition. The * 
ability of these new techniques to grow materials of alternating layers of precise 
thickness of several nanometers has totally revolutionized the field of solid state 
physics. We are nearing the age of designer materials where one grows a material to 
specifications rather than hunts for a natural material of desired properties. At present, 
commercial application of these new areas has been slow, partly because exploitation 
of current technology, particularly in software, lags device production. 
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APPENDIX 


Some Additional Band Structure Results 

As in Kane’s model applied to ll-VI materials, for purposes of calculating the band 
structure it is customary to assume a virtual crystal approximation and ignore alloy 
disorder (which later can be put in as a perturbation to discuss certain effects). When 
one solves the secular equation of the Kane model, one finds the following equation 
for the conduction, light hole and split off band 

E +(A- E^E 2 -(E g A + P 2 k 2 )E-|-AP 2 k 2 * =0 


where A is the spin-orbit splitting* 10 ) (Ai 

E g is the F 6 - r 8 band gap 
P is the momentum matrix element^ 10 ). 

With the origin is chosen to be at the top of the valence band, if A » E g and Pk, and 
including heavy holes, one can show: 


v - 2 , 2 . 

r r t Tl k 1 

E = E n + — — + — 
9 2m 2 




Conduction Band 
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2 2 

E = — Heavy Holes 


2m 


hh 
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■m 
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2 2 
Hi k 


2 2 
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2m (3 E g + 3 A) 


Split-Off Band, 


(A2b) 


(A2c) 


(A2d) 


where m is the mass of a free electron. 

Knowing the E vs. k relation, the density of states per unit volume can be calculated 
from 
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which is useful for a variety of calculations, see e.g. Eqn (3). 


Finally, in a further approximation, neglecting 2m , in (A2a) we can write a 

convenient simplification for the conduction band (where now E is referred to the 
bottom of the conduction band). 


E 



where 



(A4) 


Since the momentum effective mass is given by 
we find with Eqns. (A4) and (A5), 


(A5) 


m D 2E 
— £=1 +~ 


which shows how m p varies with energy. In many applications, only the electrons near 
the Fermi energy are important. From Eqn (A6) with E equal the Fermi energy, we see 
how m p could vary if the Fermi energy varies due to change in electron concentration. 
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FIGURE CAPTIONS 


FIGURE 1 
FIGURE 2 
FIGURE 3 
FIGURE 4 

FIGURES 
FIGURE 6 

FIGURE 7 
FIGURE 8 


The zinc-blende lattice for MCT, if the A sites are Te the B sites are 
either Hg or Cd. The unit cells have side a. 

The first Brillouin Zone of the fee lattice with several points of high * 
symmetry given the conventional labeling. 

Sketch of calculated band structure of CdTe and HgTe for higher 
energies. 

(a) Energy gap versus x for T = 300 (upper line) and T = 77 
(lower line). 

(b) Cut off wave length vs. E g . The positions of several NGS 
are listed. 

Schematic of Kane model band structure near k = 0. Results apply 
to MCT for x> 0.15. 

MCT (upper) and MZT (lower) electron mobility for the same 
energy gap. N d = 5 x 10 14 /cm 3 , N a = 2x 10 14 /cm 3 , x(MCT) = 0.193, 
x(MZT) = 0.1315. 

Schematic of different typical point defects, (a) interstitial atom - 
self or foreign, (b) Substitutional impurity, (c) Vacancy - 
absence of atom. 

Deep levels in MCT and MZT 

(a) Substitutional cations in MCT. 

(b) Substitutional cations in MZT. 

(c) Interstitials in MCT. 

(d) interstitials in MZT. 
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Schematic of Kane Energy 
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TABLE I 


Narrow Gap M-V Materials 
and 

Room Temperature Bandgap 



• E g 

InAs 

0.354 eV 

inSb 

0.180 eV 

GaSb 

0.70 eV 
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TABLE \\ 



il-V! Materials 


MCT 

HgCdTe 

Hg,. x Cd x Te 

MZT 

HgZnTe 

Hg,_ x Zn x Te 

MZS 

HgZnSe 

Hgi- x Zn x Se 
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TABLE II! 


Brief Properties of MCT 


A pseudo binary JI-VI compound with structure isomorphic to zinc blende. 

Hgi_ x Cd x Te forms a continuous range of solid solutions between the semimetal 

HgTe and the semiconductor CdTe. The bandgap is tunable from 0 to about 
1.6 eV as x carries from about 0.15 (at low temperature) to 1.0. The bandgap 
also depends on temperature, increasing linearly with temperature for a fixed 
value of x. 

Useful as an infrared detector at liquid nitrogen temperature in the wavelength 
8-12 microns which is an atmospheric window. A higher operating temperature 
than alternative materials and MCT has high detectivity, fast response, high 
sensitivity, 1C compatible and low power. 

The band structure involves mixing of unperturbed valence and conduction 
band wave function as derived by the Kane theory. They have non-parabolic 
bands which makes their analysis more difficult. 

Typical carriers have small effective mass (about 10' 2 free electron mass) 
which implies large mobility and enhances their value as IR detectors. 

At higher temperatures (well above 77K) the main electron scattering 
mechanism is the scattering by longitudinal optic modes. This scattering 
process is inelastic and it makes the calculation of electron mobility by the 
Boltzmann equation more difficult (non-iterated techniques for solving 
this equation do not work). 

The small bandgap and relatively high concentration of carriers makes it 
necessary to include screening in the calculation of the scattering of carriers by 
several interactions. 

Candidate for growth in microgravity in order to make a more perfect crystal. 



TABLE IV 


Some II-VI Compound Semiconductor Properties 


Atomic 

Number 

Hg 80 

Cd 48 

Zn 30 

Te 52 

Outer 

Electron 

5d 10 6s 2 

4d 10 5s 2 

3d 10 4s 2 

5s 2 5p 4 

Configuration 






Atomic 

Weight 

Density 

Melting 

Point 

Zn 

65.38 

7.14 g/cm 3 

693K 

Cd 

112.40 

8.65 

594K 

Hg 

200.59 

13.6 

234K 

Te 

127.60 

6.24 

723K 


a in A for unit cell 

ZnTe 6.089 

CdTe 6.481 

HgTe 6.461 
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TABLE V 


Definition and Significance of Deep and Shallow Levels 


Shallow Levels are defect levels produced by the long range Coulomb potential of 
defect. 

Deep Levels are defect levels produced by the central cell potential of defect. 


ENERGIES 


TYPICAL 

PROPERTIES 


DEEP 

SHALLOW 

May or May 
Not be Near 
Band Edge 

Near Band Edge 

Spectrum is 
Not Hydrogen 
Like 

Spectrum is 
Hydrogen Like 

Recombination 

Centers 

Suppliers of 
Carriers 

Compensators 


Electron-Hole 

Generators 
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TABLE VI 


Types of Deep Point Defects 

Substitutional Atoms 
Interstitial Atoms 
Vacancies 

Antisite (Compound Semiconductors) 
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